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ABSTRACT 

We have investigated the 2009 July impact event on Jupiter using the ZEUS -MP 2 three-dimensional hydro- 
dynamics code. We studied the impact itself and the following plume development. Eight impactors were 
considered: 0.5 km and 1 km porous (p = L760 g cm~^) and non-porous (p = 2.700 g cm~^) basalt impactors, 
and 0.5 km and 1 km porous (p = 0.600 g cm~^) and non-porous (p = 0.917 g cm~^) ice impactors. The sim- 
ulations consisted of these bolides colliding with Jupiter at an incident angle of 6> = 69° from the vertical and 
with an impact velocity of v = 61.4 km s~^ Our simulations show the development of relatively larger, faster 
plumes created after impacts involving 1 km diameter bodies. Comparing simulations of the 2009 event with 
simulations of the Shoemaker-Levy 9 events reveals a difference in plume development, with the higher in- 
cident angle of the 2009 impact leading to a shallower terminal depth and a smaller and slower plume. We 
also studied the amount of dynamical chaos present in the simulations conducted at the 2009 incident angle. 
Compared to the chaos of the SL9 simulations, where ^ ~ 45°, we find no significant difference in chaos at the 
higher 2009 incident angle. 
Keywords: comets: individual (Shoemaker-Levy 9) — hydrodynamics — methods: numerical 



1. INTRODUCTION 

Between 1994 July 16 and 22, fragments of the comet 
D/Shoemaker-Levy 9 (hereafter SL9) penetrated the Jovian 
atmosphere. This predicted impact gained great attention, 
with many E arth- and space-based telescopes aimed at this 
spectacle (seeHarrington et al.|2004| for a review of the event). 
Despite impact occurrences of this nature being characterized 
as highly unlikely, another object collided with Jupiter some- 
time between UT 9-1 1 on 2009 July 19, drastically increasing 
the expected collision rate of 0.5-1.0 km bodies with Jupiter 
([Sanchez-Lav ega et al.]|2010 i Unfortunately, the 2009 im- 
pact itself was not directly observed; however, it was analyzed 
through observations of th e impact's aftermath and was com- 
pared to th e SL9 impacts (|de Pater et al.|2010t [Fletcher et al. 
2010; Ha mmel et al.|2010[]nrton et al.|201 1[ [Sanchez-Lavega 
et al.|2010V Without a direct observation of the event, we use 



those of SL9 impact simulations. Lastly, the conclusions are 
contained in Section [5l 



2. IMPACT MODEL 



numerical simulations to seek a better understanding of the 
possible impact circumstances that could have produced this 
large atmospheric disturbance. 

We use the ZEUS -MP 2 hydrodynamics code to simulate 
the collisions of several different types of impactors, sam- 
pling the impact parameter space constrained by observations. 
From simulation results, we garner information about pos- 
sible penetration depths, plume development dynamics, and 
impact energies of potential 2009 impact scenarios. We also 
compare our results to numerical simulations of the SL9 im- 
pact event conducted by jKorycansky et aL| ( [2006| ) and jPalotai 
[et al.| p011| ). We cornpare energy deposition, penetration 
depth, and plume development between the SL9 and 2009 
simulations. 

This paper is organized as follows: Section [2| gives a de- 
scription of our numerical model. The results or the simula- 
tions are given in Section [3] Section [4| discusses the results, 
compares them to observations, and also compares them to 

parradpond@ gmail.com| 



As with [Korycansky et al.| ( |2006[ ) and [Palotai et al.[ ( [2011| ), 
we used ZEUS-MP 2 — a 3D, paral lel hydr odynamics code— - 
for impact simulations. S ee Hayes et al.[ |2Q06j fo r Z EUS- 
MP 2 code details, and see [Korycansky et al.[ ( |2002| ) and [Ko- 
[rycansky & Zahnle ( 2003[ ) for the modifications made to the 
ZEUS-MP code to include multiple materials. By extending 
both the simulation time and the spatial extent of the compu- 
tational grid in which the simulations are run, our ZEUS-MP- 
based code can model both the impact phase and the entry- 
response/blow out phase of the impac t (here, we follow the 
terminology of [Harrington et al.|2004| ). Hence, the impactor's 
entry into the atmosphere, impactor break up, and plume for- 
mation and development are simulated in each run. 

In the present simulations, we use the same coordinate sys - 
tem and a similar computational grid as [Palotai et al.| ( [2Qir[ ). 
A Cartesian coordinate system is used in the simulations: xi, 
the "along-track coordinate," is aligned with the impactor's 
initial trajectory; X3, the "cross-track coordinate," is perpen- 
dicular to the impactor's initial trajectory, and X2, the horizon- 
tal axis, is perpendicular to both xi and X3 . The local Cartesian 
coordinates for the Jovian reference frame are given by x, y, 
and z and are r elated to xi, X2, and X3 b y the following equa- 
tions , given by [Korycansky et al.[ ( [2006[ ) : 



X = X2^ 

3; = -xisin6>+X3COS^ 
Z=xicos6+X3sm0, 



(1) 



where is the angle of incidence, i.e., the angle between xi 
and the local vertical, z. The origins of both coordinate sys- 
tems coincide with the location of the 1-bar pressure level in 



the simulated Jovian atmosphere. Within and in close prox- 
imity to the impactor, the resolution is constant at 16 grid cells 
across the radius of the impactor (R16). The grid spacing then 
increases geometrically in each direction away from the R16 
area. The resolution is again held constant (4 km per grid cell) 
at the tail end of the grid, the area in which the impact plume 

develops and evolves. 

As described by |Korycansky et alT pOQ6| ), we use the Tillot- 
son equation of state (EOS). The Tillotson EOS was derived 
for cases requiring high- velocity impact calculations, can de- 
scribe the transition of shocked material into the vapor phase, 
but cannot represe nt a two-pha se region, i.e., where a liquid 
and gas co-exist ( Melosh|1989| ). The Tillotson EOS parame- 
ters used for the ba salt and ice impacto rs are the same as those 
listed in Table 1 of jKorycansky et aL[ (2QQ6). These EOS pa- 
rameters characterize the behavior of the different impactor 
materials. 



For possible 2009 bolides, Sanchez-Lavega et al. (2010) 
suggest ice impactor diameters of ^ 0.5 km up to ~ 1 km 
based on comparisons to SL9 models a nd ablation rate consid- 
erations at higher angles of incidence. |Hammel et al.| ( [2QT0| ) 
suggest ice impactor diameters of 500-700 m based on sim- 
ilarities of the 2009 impact site to the E and R impact sites 
of SL9. There is observational evidence suggesting the pos 



sibility that the 2009 impactor was asteroidal in origin, rather 
than cometary dFletcher et al.|2010t|Hammel et al.|2010t[Qr^ 



|ton et al.||20il ). From therm al heating and mass transport 



estimates, |Orton et al.| ( |2011| ) suggest diameters of 200-500 
meters for basalt impactors of density 2.5 g cm~^ (an im- 
pactor mass range of ~ 1.05 x 10^^ grams to ~1.64 x 10^"^ 
grams). The following eight impact cases were run to ~ 30 
seconds after impact: 0.5 km and 1 km porous (p = 1.760 g 
cm~^) and non-porous (p = 2.700 g cm~^) basalt impactors, 
and 0.5 km and 1 km porous (p = 0.600 g cm~^) and non- 
porous (p = 0.917 g cm~^) ice impactors. We model 0.5 km 
bodies in order to sample relatively smaller impactors that 
satisfy size estimates for both ice and basalt bolides, and we 
model 1 km diameter impactors for easy comparison to pre- 
vious SL9 models. An incident angle of 69° from the vertical 
and a n impact latitude of 55°. 10 S were used in the simula- 
tions ( jSanchez-Lavega et al.|2010| ). The gravitational acceler- 
ation at this latitude, including the J2 and centrifugal terms, is 
2582 cm s~^. An impact velocity of v = 61.4 km s~^ was used 
for the purposes of comparison to previous S L9 simulations 
( [Korycansky et al.||2006[ [Palotai et al.||201 1) . Several runs 
were also conducted to test the degree of dynamical chaos, 
the sensitivity of results to initial conditions, present in simu- 
lations with 2009 impact parameters ( Korycansky et al. 2006) . 

3. RESULTS 

Impact parameters for all the majo r cases run in t he pre sent 
paper, plus the parameters used by Palotai et al. ( 2011| ) for 
the SL9 case, are given in Table [l| and a summary of sim- 
ulation results is given in Table [2p A case label is given to 
each simulation. Figure [T] gives a series of snapshots of the 
1 km non-porous ice impactor' s simulated decent into the Jo- 
vian atmosphere. This figure shows the typical progression 
of our simulations. The first row of panels, ^ = 4.00 seconds 
after impact, shows the impa ctor traveling before it has be- 
gun to fall apart. As shown by JFalotai et al.| ( [20TT| ), all of the 
impactor material is contained within a narrow trail follow- 
ing the body, constrained by the shock system on the trailing 
edge of the impactor. In the second row, ^ = 8.00 seconds, the 







Table 1 










Impact Parameters 






Case 


Material 


Density 


Diam. 


Angle 


Latitude 


label 




(g/cm^) 


(km) 






I05p 


Ice 


0.600 


0.5 


69° 


55°.10S 


I05n 


Ice 


0.917 


0.5 


69° 


55°.10S 


B05p 


Basalt 


1.760 


0.5 


69° 


55°.10S 


B05n 


Basalt 


2.700 


0.5 


69° 


55°.10S 


IlOp 


Ice 


0.600 


1.0 


69° 


55°.10S 


IlOn 


Ice 


0.917 


1.0 


69° 


55°.10S 


BlOp 


Basalt 


1.760 


1.0 


69° 


55°.10S 


BlOn 


Basalt 


2.700 


1.0 


69° 


55°.10S 


^SL9p 


Ice 


0.600 


1.0 


43°. 09 


44°.02S 


SL9n 


Ice 


0.917 


1.0 


43°. 09 


44°.02S 


^ SL9 parameters used 


in|Palotai et a 


.|l201l|. 







body has begun significant breakup and decelerates. The im- 
pactor material near and around the body begins to spread out 
quickly, and the shock system becomes turbulent. Material 
left behind by the bolide in the back half of the grid, within 
the high- temperature, low-density region that will form the 
plume, spreads out less rapidly. By ^ = 12.00 seconds, the 
impactor has become very incoherent and has almost reached 
its terminal depth. Impactor material and Jovian air in the 
back region begin their rapid ascent in the atmosphere, signi- 
fying the start of plume genesis. The plume is indicated by 
the growing red region. The majority of the impactor mate- 
rial travels to its terminal depth between t = 16.00 and 24.00 
seconds, and the plume can be seen speeding up and rising 
straight up in the X3 direction. As the impactor travels deeper, 
its vapor diffuses into the Jovian atmosphere, and the impactor 
fraction decreases with time. The plume increases its blowout 
speed and rises well above the initial impact trajectory of X3 
= 0. The impactor material in this region also continues to 
dissipate as the plume travels in the atmosphere. In general, 
all the impact simulations presented in this paper proceed as 
described above and resemble that pictured in Figure [T] Note- 
worthy differences in the depths reached by the impactors and 
the development of each plume do exist, however, and are ex- 
plained in the following subsections. 

3.1. Terminal Depth and Energy Deposition 

Figure [^ shows the kinetic energy deposition curves for the 
porous and non-porous 2009 impact cases. Also included for 
comparison on each plot is the energy deposition curve of the 
appropriate SL9 impact case. These plots, similar to those of 
[Korycansky et a l. ( 2006 ), show the amount of energy an im- 
pactor releases to the surrounding atmosphere per unit of alti- 
tude traveled. The location of the initial and sudden increase 
of an energy deposition curve is an indicator of the height at 
which rapid and extensive structure loss of the impactor be- 
gins. Maximum energy deposition occurs when the nucleus of 
the impactor loses all coherency, and the impactor explodes. 

The altitude at which the energy deposition returns to a 
value of zero is the terminal depth of the impactor. Each of 
the simulated impactors begins significant breakup at different 
altitudes and reach varying terminal depths. These terminal 
depths are listed in Table [2] and are given in the local vertical 
coordinate, z. A positive z represents an altitude above the 1- 
bar level in the Jovian atmosphere, and a negative z represents 
an altitude below the 1-bar level. The terminal depths are 
consistent with energy considerations: larger, more dense im- 
pactors will penetrate the deepest, whereas smaller, less dense 
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Figure 1. Progression of our impact simulations. Left: Along-track velocity values [km s~^ ], i.e., the velocity in the x\ direction. A red color indicates an upward 
velocity, material moving higher into the Jovian atmosphere; a blue color indicates a downward velocity, material moving deeper into the atmosphere. The red 
areas are indicators of the rising plumes. Right: the fraction of impactor material compared to Jovian air. The along-track coordinates are given at the bottom of 
the plots. The IlOn case is shown here. The first panel shows an arrow indicating the direction of up in local Jovian coordinates. Note: the x\ and X3 length scales 
are not the same. 

moves deeper down the length of the plume. The ejection an- 
gles of the plumes are similar across all the 2009 cases, and 
within the simulated ~ 30 seconds, the plumes attain ejection 
angles of 60° -70° from the vertical. 



impactors will reach shallower depths. The total energies of 
the impactors, obtained by integrating dE/dz, are given in Ta- 
bleEl 

3.2. Plume Development and ''Pinch-off" Regions 

Figure [3] shows the plume structure at comparable times 
of plume evolution for each porous 2009 impact simulation. 
Figure |4] shows the same for the non-porous 2009 impact 
simulations. Within the simulated ~ 30 seconds, higher re- 
sultant ejection speeds are reached in the 1 km diameter 
cases. Markedly smaller plumes and generally lower ejections 
speeds are seen in the 0.5 km impactor cases. 

Figure |5] shows the velocity distribution of a rising plume 
and gives the fraction of impactor material contained within. 
All of the plume ejecta are above the initial impact trajectory. 
This is different from the SL9 simulations, where the plume 
expands to a greater diameter in the simulated ti me and still 
crosses the initial impact path ( [Palotai et al.|2011| ). The distri- 
bution of impactor material w ithin the plume issim ilar to the 
SL9 simulation conducted by [Palotai et al. ( 2011 ), however. 
A maximum of ~ 10% impactor material is located near the 
top of the ejecting plume, and this fraction decreases as one 



Just as in the SL9 simulations conducted by Palotai et al. 
( |2011| ), a pinch-off region appears in each of the 2009 cases. 
These regions are shown in Figures [3] and |4] Above this 
pinch-off level, heated Jovian atmosphere containing a small 
fraction of impactor material rises and expands as a plume 
(shades of red in Figures |3] and |4l; below this level, a major- 
ity of the impactor material continues downward (shades of 
blue in Figures |3] and H]) and will later rise more slowly and 
independently ofthe ejecting plume as a bubble-like region of 
impactor material. The approximate location of this pinch-off 
region for each case is listed in Table [^ 

The altitudes at which the pinch-off regions occur appear 
to correspond with the altitudes at which the impactor b egins 
rapid and extensive structure loss ( [Palotai et al.||201l) . We 
are currently working to better characterize and garner ad- 
ditional details about the relationship between the pinch-off 
region and impactor break-up. 

4. DISCUSSION 



Table 2 

Simulation Results 



Case 
label 



Terminal Depth^ 
(km) 



Pressure at 
Terminal Depth 

(bar) 



Total Energy 
(erg) 



Maximum Plume Velocity 
(km/s) 



Pinch-off Location^ 
(km) 



Pressure at 

Pinch-off Location 

(bar) 



I05p 


-24 


2.44 


7.4 X 1026 


7.5 


I05n 


-36 


3.52 


1.1 X 1027 


9.1 


B05p 


-54 


5.67 


2.2 X 1027 


8.0 


B05n 


-55 


5.83 


3.3 X 1027 


8.3 


IlOp 


-61 


6.72 


5.9 X 1027 


10.3 


IlOn 


-80 


10.2 


9.1 X 1027 


12.5 


BlOp 


-102 


15.8 


1.7 X 1028 


11.0 


BlOn 


-114 


19.7 


2.6 X 1028 


11.2 


bSL9p 


-124 


23.5 


5.9 X 1027 


17.2 


SL9n 


-150 


35.0 


9.1 X 1027 


16.0 



21 


0.365 


19 


0.408 


-1 


1.06 


-2 


1.10 


4 


0.891 


-4 


1.19 


-18 


2.00 


-23 


2.37 


-33 


3.23 


-50 


5.15 



^ Given in z, altitude relative to the 1-bar pressur e level. 

Results from simulation using parameters from|Palotai et al.|J201 
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Figure 2. Kinetic energy deposition curves for the porous impact simulations 
(a) and the non-porous impact simulation (b). This plot demonstrates the 
increase of both terminal depth and total energy deposited as the size and 
density of the impactors grow. 



4. 1 . Observations and Simulations 



Ob s ervations conducted by Jde Pater et al.|(|201Q|), [Fletcher 
eraLl([2QTQl),|Hammel et al.l ( |201U| ), |Orton et al.| ( |201 1D, and 



Sanchez-Lave ga et al.| (|2010) place several constraints on the 
impact parameters and the dynamics of the 2009 impact. In 



this section, we will examine these constraints and compare 
them to the results of the present simulations. 

Two interesting constraints include the extent of debris dis- 
tribution and temperature perturbations in the Jovian atmo- 
sphere as a result of the impact. Debris deposition was con- 
strained rou ghly between the 10-mbar and 300-mbar level s 
dde Pater et al. 2010; Ha mmel et al.|2010l[Qrton et al.|2011| ). 
Thermal perturbations reached higher up in the atmosphere 
than the debris, up to about the 0.1-20-mbar levels, and were 
present slightly deeper into the atmosphere than the de bris, 
down to the 400-600 mbar levels (de Pater et al. 2010; Ofton] 



et al.|20lT] ). A lack of excess methane emission in the upper 



stratosphere is also an ind icator of heating lim ited to pres- 
sures greater than 10 mbar ( Fletcher et al.l2010 ). This is very 
different from the the SL9 impacts, in which plumes reached 
thousands of kilometers above the 100-mbar level in Jupiter's 



atmosphere ( [Hammel et al.|1995} Jessup et a l. 2000). 

Enhanced levels of ammonia in the lower stratosphe re also 
place boundaries on the penetration depth of the plume. [Orton 



et al.[([20lT|) found that a vertical profile of ammonia peaking 



in the 20-36-mbar region is required to reproduce the spectral 
shape of the 2009 NH3 emission at the impact location. This 
enhanced presence of ammonia in the Jovian stratosphere im- 
plies the 2009 impact wake probably reached down to the 600- 
700 mbar levels, retrieving tropospheric ammo nia and trans- 
porting it to the stratosp here as the plume rose ( |de Pater et al. 



2010t [Orton et al.|2011| ). However, it is likely that the down- 
ward et that would form the plume did not penetrate mucl: 



jet that would form the plume did not penetrate much 
farther than the 700-mbar level since the NH3 gas was con- 
tained close to th e center component of the impact streak ( [de 
Pater et al.|2010 ), and the jet probably did not make it to trie 



Jovian water cloud level since shock chemisty of the 2009 
event favored prod uction of ethane an d other hydrocarbons 
over CO and H2O ( (Fletcher et al.|2010| ). 

Table [2] gives the pressure levels at which the pinch-off re- 
gions occur for each of the conducted simulations. The pinch- 
off region is a measure of the deepest penetration of the fire- 
ball that forms the plume and is the lower limit from which 
material may be dredged up from the Jovian depths. In gen- 
eral, the plume jet is small and weak at the pinch-off region, 
but grows in speed as one moves up the plume channel. Cases 
I05p and I05n have pinch-off regions located at ~ 3 65 -mbar 
and ~4 10-mbar, respectively. Ice impactors of these sizes, 
densities, and impact velocity result in plumes that do not 
penetrate down to the stratospheric ammonia reservoir. Cases 
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Figure 3. Plume structures of the porous 2009 impact cases. Along-track velocity is shown. The same color bar used for the left panels of Figure^ also 
applies here. Top-left: 0.5 km, porous ice; Top-right: 0.5 km, porous basalt; Bottom-left: 1 km, porous ice; Bottom-right: 1 km, porous basalt. The along-track 
coordinates are given at the bottom of the plots. The top of the plots indicates lines of constant height, z [km], in Jupiter's atmosphere, z = km represents the 1 
bar pressure level. Note: the xi and X3 length scales are not the same. 



_ 20r 
^ 10- 



Ice; p = 0.917 g/cm^ 
20 40 



60 



80 



-20 



Basalt; p = 2.700 g/cm^ 
20 



40 



60 



-10 



-20 



Pinch-Off Region 



Impactor Entry Path 

I05n; Time: 20.00s 



Pinch-Off Region 



B05n; Time: 20.20s 



100 



-50 



50 



100 



150 



-150 



-100 



-50 



50 



100 



20r 



-60 



-40 



-20 



20 



40 



60 



80 



-80 



-60 



-40 



-20 



20 



40 



\ V 



\ \ 



■X — \ — \ — X 



10- 



-10 



Pinch-Off Region 



llOn;Time: 21.20s . . 



Pinch-Off Region 



BlOn; Time: 20.60s 



-20 



-200 



-100 100 

x^, Along-track Coordinate (km) 



200 



-300 



-200 -100 

x^, Along-track Coordinate (km) 



IOC 



Figure 4. Plume structures of the non-porous 2009 impact cases. Along-track velocity is shown. The same color bar used for the left panels of Figure^ also 
applies here. Top-left: 0.5 km, non-porous ice; Top-right: 0.5 km, non-porous basalt; Bottom-left: 1 km, non-porous ice; Bottom-right: 1 km, non-porous basalt. 
The along-track coordinates are given at the bottom of the plots. The top of the plots indicates lines of constant height, z [km], in Jupiter's atmosphere, z = km 
represents the 1 bar pressure level. Note: the xi and X3 length scales are not the same. 



B05p and B05n result in relatively small plumes whose pinch- 
off regions penetrate down to the ~1.10-bar level, past the top 
levels of Jupiter's NH3 clouds, and speeds in the lower plume 
jet around the 700-mbar level reach ~1.8 km s~^ for B05p 
and '^1.0 km s~^ for B05n. Cases IlOp and IlOn have pinch- 
off regions located at ~890-mbar and ~1.2-bar, respectively. 
For these two 1 km ice cases, speeds in the lower plume jet at 
~700-mbar reach about 1.8 km s~^ for IlOp and 2.8 km s~^ for 
II On. Cases B lOp and B lOn result in larger plumes that pene- 
trate well past 700-mbar, down to pressures around the 2 and 
3 bar levels, respectively, but stay above the H2O cloud tops 
around 6 bar. At the ~700-mbar level, speeds in the lower 
plume jet reach about 4.6 km s~^ for BlOp and about 5.6 km 
s-iforBlOn. 

Cases I05p and I05n do not satisfy the observational con- 
straints, so it is likely that bodies with these impact character- 
istics did not cause the atmospheric response seen on Jupiter 



in July 2009, setting a lower limit on the possible size, density, 
and impact velocity of the bolide. The plumes of cases B05p, 
B05n, IlOp, IlOn, BlOp, and BlOn penetrate deep enough 
to reach the ammonia clouds in the Jovian atmosphere, and 
so these impactors remain possible 2009 candidates. Further 
constraining the possibilities of the 2009 impact will require 
continued exploration of the possible impact parameter space, 
detailed ammonia transport calculations, and extending simu- 
la tions both spatially and temporally. 



Orton et al 



( |2011 ) give 7 ± 2 X 10^^ erg as a lower esti- 



mate for the energy of the impact. This is a lower estimate 
because it does not take into account the 4 days of cooling 
that passed between the impact and the observations nor the 
large amounts of energy lost to other dynamical processes in 
the impact, such as plume formation and the transport of at- 
mospheric waves. All of the present simulations satisfy this 
lower limit. 
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Figure 5. Velocity vector field of the impact plume for case IlOp, 27.0 seconds after impact. Color gives the fraction of impactor material present in the rising 
plume. The straight, dark line indicates the initial trajectory of the impactor. 
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Figure 6. Comparison of the atmospheric density distribution during the aftermath of two impact events, one at the SL9 impact angle (left) and one at the 2009 
impact angle (right). The undisturbed Jovian atmospheric density profile is plotted in the background. Contours give the log(p) values for each impact and 
Jupiter's atmosphere. The local vertical is given on the left axis; the corresponding pressure values are given on the right axis. For clarity, the SL9 case is shifted 
to the left in the _y coordinate by 400 km relative to the 2009 case. Both impactors are identical 1 km non-porous ice bodies impacting at 61.4 km s~^ (cases SL9p 
and IlOp). The straight, dark lines indicate the initial trajectories of the impactors and the dotted lines indicate the boundaries of the simulations' computational 
grids. 

both impactors are identical and share general characteristics, 
the terminal depths and plume development are strikingly dif- 
ferent at the differing incident angles. 

In this particular comparison, the terminal depth of the SL9 
impactor is ~ 60 km deeper than the 2009 impactor' s terminal 
depth. This is mainly due to the larger incident angle of the 
2009 impact: for a given distance traveled along the path of 
each impactor (in xi), displacement in height z will be smaller 
at higher incident angles. The pinch-off regions differ by a lit- 
tle more than a scale height, with the SL9 pinch-off reaching 
to ~ 30 km below the 1-bar level and the 2009 pinch-off oc- 
curring ~ 4 km above the 1-bar level. There is also a distinct 
difference between the developing plumes. For the SL9 case, 
the plume begins to expand significantly as it rises in the at- 
mosphere. The SL9 plume reaches a diameter of ~35 km by 
30 seconds after impact. Such rapid expansion is not seen in 
the 2009 plume, however. During the formation of the 2009 
plume, it rises above the impact path while undergoing rela- 



4.2. 2009 vs. SL9 Impacts 

Figure [6| shows a snapshot of the aftermath of both an SL9 
impact and a 2009 impact. The figure gives the density distri- 
bution for both cases, plotted against the undisturbed Jovian 
density distribution. Case SL9p and case IlOp are used in 
the figure, i.e., the impacts in this figure are identical, except 
for the incident angle and the latitude at which the impactor 
strikes Jupiter. Both cases exhibit features characteristic of all 
our present impact simulations. At the terminal depths of the 
impactors, a relatively-low-density region containing most of 
the impactor material exists. This is the bubble-li ke region be 



low th e pinch-off level seen in the simulations of Palotai et al. 
( |2011| ). Shocks can be seen propagating away from the impact 
path, and at the boundaries of plume formation, two shock 
waves can be seen. The lowest-density region marks the loca- 
tion of the plume. In the upper atmosphere, the plume gains 
speed, rises, and expands; in the lower atmosphere, the plume 
becomes slender, tapering to the pinch-off region. Though 
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Figure 7. Dynamical chaos analysis. Kinetic energy deposition curves at the 
2009 incident angle (a) and the SL9 incident angle (b). The black line gives 
the energy deposition for the nominal case of a 1 km porous ice impactor 
traveling at v = 61 .4 km s~^ . The bottom plot is Figure 7 (d) from Korycansky 
|itdH2006t . ' 

tively little expansion. By 30 seconds after impact, the 2009 
plume only reaches a diameter of ~15 km. The difference 
in plume sizes can be seen in Figure [6] which shows the SL9 
plume about twice the size of the 2009plume. Table [2]reveals 
that the SL9 plume contains air and entrained impactor ma- 
terial traveling at speeds about 7 km s~^ greater compared to 
the 2009 plume. The larger plume and heightened speeds in 
the SL9 case relative to that of 2009 implies a dependance of 
plume size and maximum plume height on incident angle. 

4.3. Dynamical Chaos 

In order to observe the sensitivity of the present results 
to initial conditions, we conducted several 1 km porous ice 
impactor simulations with slightly different initial conditions 
compared to a nominal case. These differences in initial con- 
ditions include a change in the impact velocity (Av) by 0.1% 
and a shift in the initial position of the impactors (in each di- 
rection xi , X2, and X3) by a half a grid cell, ~ 15 m in this case. 
The variation in terminal depths and peak energy deposition 
locations is a measure of the dynamical chaos present in the 
simulations. Korycans ky et al., ( ,2006j ) conducted this same 
analysis and found that among the four SL9 impact cases they 
simulated, a 1 km porous ice impactor was subject to the high- 
est degree of dynamical chaos. Only slightly changing initial 



conditions for this case resulted in relatively large variations 
in the terminal depths and peak energy deposition locations. 

Figure [7] (a) gives our chaos analysis of the 2009 impact 
for the porous ice case. Comparing Figure [7] to Korycansky 
et al.'s (2006) Figure 7 (d), seen here as our Figure [7] (b), we 
do not find the dynamical chaos to be significantly aifferent 
at an incident angle of 6> = 69°. The standard deviation of the 
terminal depths at the 2009 impact angle is 5'/o9 ~ 13.2 km in 
the along-track coordinate, and the range in terminal depths 
is about 44 km. At the 6> ?^ 45° SL9 angle of incidence, the 
standard deviation of the terminal depths is Ssu^ ^20.8 km in 
the along-track coordinate, and the range in terminal depths is 
about 53 km. 

We conducted a simple F-Test to compare the variation in 
terminal depths between the 2009 and SL9 cases. Our test F 
statistic, F - S\^^lS\^ - 2.48, was compared to the critical F 
statistic at the 5% level, F = 4.39. Though visually the termi- 
nal depths for the 2009 porous ice case seem less scattered, 
and the range and standard deviation of these terminal depths 
is less at the 2009 incident angle, we cannot say that the stan- 
dard deviation of the terminal depths at the SL9 incident angle 
is significantly larger than at the 2009 incident angle because 
our test F statistic was not greater than the critical value. 

5. CONCLUSIONS 

We have presented several possibilities for the impact and 
immediate aftermath of the 2009 collision into the Jovian at- 
mosphere. The results of the present simulations provide in- 
sight into the impact event and also provide information about 
the variation in atmospheric response due to changes in im- 
pact parameters. At the estimated 2009 incident angle of 
6/ = 69°, we see several differences between plumes gener- 
ated from 0.5 km impactors and 1 km impactors. Within the 
simulated r^ 30 seconds, the 0.5 km impactor events produce 
relatively smaller and slower plumes while the 1 km impactor 
events produce relatively larger and faster plumes. 

The penetration depths of the impactors and the pinch-off 
regions are associated with the nature of the impactors: at a 
given incident angle, the larger the impactor and the heavier 
the impactor material, the deeper the locations of the terminal 
depth and the pinch-off region. Dynamical chaos present at 
the 2009 incident angle for the 1 km porous ice impactor did 
not prove to be signific antly less than that of the most chaotic 
impactor case given by 'K orycansky et aL (20061. 

Comparing the aftermaths of an SL9 impact and the 2009 
event reveals several differences that may have consequences 
for the observable manifestation of Jupiter's atmospheric re- 
sponse. The impact plume produced at the SL9 incident angle 
is significantly larger and faster than that of the impact plume 
produced at the 2009 incident angle. 

Given observations of thermal perturbations, debris depo- 
sition, and ammonia transport in the Jovian atmosphere after 
the 2009 event, constraints have been placed on the possible 
outcomes of the im pact, including plume speeds and pinch- 
off region locations (|de Pater et al.|20Tm [Fletcher et al.|2 010 



Hammel et 



aT][2 
m th 



2010t lOrton et al.|[20l7| . Of the eight cases 
considered in the present paper, the 0.5 km ice cases cannot 
explain the atmospheric disturbance observed after the 2009 
impact; these impactors' plume jets do not penetr ate deep 
enough to explain stratospher ic NH3 observations (de Pater 
[eTaLllMTOl lOrton et al.][2QTT] ). We thus set a lower hmit on 
the size and density of the 2009 impactor. All 1 km impactor 
plumes and the 0.5 km basalt impactor plumes reach down 
to the ammonia ice cloud level in Jupiter's troposphere. To 



better constrain the impact characteristics of the 2009 event, 
more simulations are required spanning more of the possi- 
ble parameters proposed by observations. Ammonia transport 
from the upper troposphere to the stratosphere by means of 
the rising plume jets must also be looked at in more detail. 
Through this analysis, we will narrow the range of possible 
impactor circumstances that produced the atmospheric distur- 
bance on Jupiter in 2009. 
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